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. We analyze the reaction dynamics of central Pb+Pb collisions at 160 GeV/nucleon. First 



we estimate the energy density e pile-up at mid-rapidity and calculate its excitation function: 



e is decomposed into hadronic and partonic contributions. A detailed analysis of the collision 
dynamics in the framework of a microscopic transport model shows the importance of partonic 
degrees of freedom and rescattering of leading (di) quarks in the early phase of the reaction for 
Ei a b > 30 GeV/nucleon. The energy density reaches up to 4 GeV/fm 3 , 95% of which are contained 
in partonic degrees of freedom. It is shown that cells of hadronic matter, after t ~ 2R/jv cm , 
can be viewed as nearly chemically equilibrated. This matter never exceeds energy densities of 
~ 0.4 GeV/fm 3 , i.e. a density above which the notion of separated hadrons loses its meaning. The 
final reaction stage is analyzed in terms of hadron ratios, freeze-out distributions and a source 



analysis for final state pions. 
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The study of relativistic heavy ion collisions offers the unique opportunity to study hot and dense 
QCD matter under conditions which are thought to have existed in the early stages of our universe. 
However, only the hadronic final state of the heavy ion collision is accessible via experiment, or - in the 
case of leptonic probes - the time integral of the emission over the entire reaction history. Microscopic 
transport models offer the unique opportunity to link this final state information to the experimentally 
inaccessible early and intermediate reaction stages. In this paper we analyze the reaction dynamics of 
central Pb+Pb collisions at CERN/SPS energies. We focus specifically on the time evolution of energy 
density and its interpretation in terms of hadronic and partonic degrees of freedom. We then discuss 
the possible formation of a thermally and chemically equilibrated state in the central reaction zone 
and finally investigate the late reaction stages with a decomposition of freeze-out radii and sources for 
individual hadron species. 

The determination of energy densities in ultra-relativistic heavy-ion collisions is crucial for any 
discussion involving a possible deconfinement phase transition to a QGP [p], 0, [|, |5], || 0] . Estimates 
for the energy density during the hot and dense early reaction stage have been given by a large variety 
of different models g, 1, g g 1 0. 

It has been questioned whether hadronic transport models are still valid at CERN/SPS energies: 
the energy density estimates obtained in these frameworks are claimed to be well above the critical 
energy density estimates for a deconfinement phase transition given by Lattice Gauge Theory Jll], 
PI [T^l . Hadronic transport models, however, contain implicit partonic degrees of freedom - particle 
production at high energies is e.g. modeled via the production and fragmentation of strings |L4], [15], |X6f . 

In the UrQMD model used below, the leading hadrons of the fragmenting string contain the 
valence-quarks of the original excited hadron. These leading hadrons are allowed - in the model - to 
interact even during their formation time, with a reduced cross section, thus accounting for the original 
valence quarks contained in that hadron. Those leading hadrons represent a simplified picture of the 
leading (di) quarks of the fragmenting string. Newly to-be-produced hadrons which do not contain 
string valence quarks do in the present model not interact during their formation time - however, 
they contribute to the energy density of the system. A proper treatment of the partonic degrees 
of freedom during the formation time ought to include soft and hard parton scattering |28| and the 
explicit time-dependence of the color interaction between the expanding quantum wave-packets p5| : 
However, such an improved treatment of the internal hadron dynamics has not been implemented for 
light quarks into the present model. Therefore, in the following analysis all contributions stemming 
from hadrons within their formation time are termed "partonic". All contributions stemming from 
fully formed hadrons are termed "hadronic" . The main focus of this paper is on the partitioning and 
the time evolution of the energy density and the collision dynamics of the early, intermediate, and late 
reaction stage at energies E lab = 10 — 200 GeV/nucleon. 



The UrQMD model [ IT | is based on analogous principles as (Relativistic) Quantum Molecular 
Dynamics |T8], [H5], [2D], . Hadrons are represented by Gaussians in phase space and are propagated 

according to Hamilton's equation of motion. The collision term of the UrQMD model treats 55 different 
isospin (T) degenerate baryon (B) species (including nucleon-, delta- and hyperon- resonances with 
masses up to 2 GeV) and 32 different T-degenerate meson (M) species, including (strange) resonances 
as well as their corresponding anti-particles, i.e. full baryon-antibaryon symmetry is included. Isospin 
is treated explicitly. For hadronic excitations with masses m > 2 GeV (B) and > 1.5 GeV (M) a 
string model is used. Particles produced in the string fragmentation are assigned a formation time. 
This time Tf physically consists of a quantal time tq, i.e. before the partons are produced, tq ~ 1/m, 




Figure 1: Time evolution of the energy density e in central Pb+Pb reactions at 160 GeV/nucleon. 
e has been decomposed into "partonic" and "hadronic" contributions and only particles around mid- 
rapidity have been taken into account. The early and intermediate reaction stages are dominated by 
the "partonic" contribution. 

and a quantum diffusion time, 7~d, during which the partons evolve in the medium to build up their 
internal asymptotic wave-functions to form the hadron. tq and tjj differ for different parton and 
hadron species. For our present purpose, we - for the sake of simplicity - just collect all partons, 
formed and unformed, as one species. For a detailed overview of the elementary cross sections and 
string excitation scheme included in the UrQMD model, see ref. |L7] . 



The partitioning of the distinct constituents can be inspected in Figure [I] which shows the time- 
evolution of the energy density for central Pb+Pb collisions at 160 GeV/nucleon. The nuclei are 
initialized such that they touch a t = fm/c. The energy density is partitioned into the above defined 
"hadronic" contribution, from fully formed hadrons, and the "partonic" contribution, from partons, 
constituent quarks and diquarks within the hadron formation time. Nearly all incident baryons are 
rapidly excited into strings. Subsequently, "partonic" energy density builds up, reaching values of 4 
GeV/fm 3 around midrapidity, Ay — 1 (at t fa 1 fm/c). In the course of the reaction hadrons are formed 
which increases in the "hadronic" energy density, accompanied by a nearly exponential decrease in the 
"partonic" energy density. 

These energy densities are calculated as follows: In the UrQMD model hadrons are represented 
by Gaussian wave packets. The width of the Gaussians o = 1.04 fm and their normalization are chosen 
such that a calculation of the baryon density in the initial nuclei yields ground state nuclear matter 
density. The (energy-) densities in the central reaction zone are obtained by summing analytically 
over all Gaussian hadrons around mid-rapidity (y c .m. ± 1) at the locations of these hadrons and then 
averaging over these energy densities. This summation over Gaussians yields a smooth estimate for 
baryon- and energy-densities, as compared to counting hadrons in a test volume. The rapidity cut 
insures that only those particles are taken into account which have interacted. Thus, the free streaming 
"spectator" matter is discarded. 

The absolute value of the energy density, however, may depend on the rapidity cut: Without 
rapidity cut the energy densities during the early reaction stage (t ~ 1 fm/c) can be as high as 20 
GeV/fm 3 . Even higher values in e can be obtained by choosing the geometric center of the collision 




Figure 2: Top: Time evolution of the multiplicity of hadrons and partonic constituents, divided 
into baryonic and mesonic contributions. Bottom: Collision rates for baryon-baryon (BB) and meson- 
meson (MM) collisions. The rates have been decomposed into interactions involving formed hadrons 
and those involving partonic constituents. 

for the sum over the Gaussians instead of averaging over the energy densities at the locations of the 
hadrons. The energy density at a single point may not be physically meaningful and therefore the 
latter method is favorable. 

The time evolution of partonic constituents and hadrons is shown in the upper frame of figure |2|. 
The first 5 fm/c of the reaction are dominated by the partonic constituents. The long-dashed and the 
dotted curves show the number of baryons and mesons contained in those constituents. In the case of 
leading-particles these can be interpreted as constituent (di) quarks or, for freshly born partons with 
small cross sections, as excitation modes of the color field. 

The lower frame of figure |2| shows the time evolution of the number of baryon-baryon (BB) and 
meson-meson (MM) collisions, both for "hadronic" and "partonic" interactions. "Partonic" interac- 
tions denote interactions of leading (di) quarks either among themselves or with fully formed hadrons. 
The early reaction stages, especially the MM case, is clearly dominated by those "partonic" interac- 
tions. This number increases further if the scattering of the newly formed partons is included. Thus 
"partonic" degrees of freedom significantly contribute both, to the energy density, as well as to the 
collision dynamics in the first 5 fm/c. 

It should be noted that the "partonic" collision rates can increase with the partonic cross section 
during formation time: In this analysis all interactions during formation time have been considered 




Figure 3: Top: excitation function of the maximum total energy density mid-rapidity hadrons expe- 
rience. Bottom: excitation function of the maximum "partonic" fraction of energy density Already 
at a beam energy of 40 GeV/nucleon more than 90% of the energy density is contained in partonic 
degrees of freedom at one time during the collision. 

purely "partonic". Other scenarios, however, include a "hadronic" contribution to the cross section 



which increases continuously during td and reaches its full hadronic value at the end of td IP5 |. 

Do "partonic" degrees of freedom play any role at 10 GeV/nucleon, i.e. at the AGS? The upper 
frame of figure |3] shows the maximum total energy density obtained in central collisions of heavy nuclei 
as a function of incident beam energy, starting from 2 GeV/nucleon and going up to 200 GeV/nucleon. 
The energy density is obtained by the same method as used figure p]. However, here "partonic" and 
"hadronic" contributions have been summed, e increases monotonously with the beam energy, reaching 
values as high as 4 GeV/fm 3 for SPS energies, which would seem unreasonably high, if a purely hadronic 
scenario were used. 

The lower frame of figure |3] shows the maximum fraction of the energy density which is contained 
in "partonic" degrees of freedom. Even at AGS, energies already more than half of the energy density 
is due to such "partonic" degrees of freedom, even though these do not yet dominate the "hadronic" 
contributions. At 40 GeV/nucleon, the maximum of the fraction of "partonic" energy density is already 
> 90% of the total e. 

The monotonous increase of the energy density excitation function does not imply that the 
excitation function of the space-time volume of high baryon density shows the same behavior. At 
AGS energies, Ei^ ~ 10 GeV/nucleon, baryons still dominate the composition of the hadronic mat- 
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Figure 4: Temperature versus baryonic chemical potential predicted by the statistical model for the 
time evolution of the hadronic characteristics {e, ps, Ps} obtained within UrQMD in the central cell 
(2x2x1 fm 3 cell) of central A+A collisions at 10.7, 40 and 160 GeV/nucleon. the open symbols 
show the non-equilibrium stage of the reaction, whereas the full symbols denote the phase of local 
thermal equilibrium in the central cell. The two solid lines correspond to the phase-boundary between 
a confined and deconfmed phase calculated for two different bag constants, B 1 ^ = 227 and 302 MeV 
corresponding to T c =150 and 200 MeV at ps = 0. 

ter, whereas at CERN/SPS energies, 200 GeV/nucleon, mesons constitute the largest fraction of the 
hadronic matter. The maximum space-time volume of dense baryonic matter can be reached at beam 
energies around 40 GeV/nucleon. A detailed analysis of that regime, also with respect to experimental 
signatures, is presently underway . 

The time evolution of temperature T versus baryonic chemical potential pb for Au+Au reactions 
at 10.6 GeV/nucleon and Pb+Pb reactions at 40 and 160 GeV/nucleon, respectively, is plotted in 
figure |j. The thermodynamic quantities T and ps have been extracted by fitting a statistical model 
to the quantities energy-, baryon- and strangeness density {e, pb, ps} obtained from UrQMD in the 
central cell (2x2x1 fm 3 cell) of the heavy-ion reaction [p^| . Here, the densities are calculated by 
summing over all relevant partonic and hadronic degrees of freedom in the cell. We see that the average 
Pb in the reaction drops drastically with the initial collision energy, while the maximal temperature is 
growing and practically reaches the upper phase transition boundary with the critical temperature of 
T c =200 MeV, as calculated with the MIT bag model (details of the used bag model can be found in |24j ). 
However, during the early reaction stages matter in the central reaction cell is neither fully hadronic, 
nor thermally and chemically equilibrated. A detailed analysis of velocity distributions and particle 
spectra in the central cell p3| reveals that at approximately t = 2 fm/c the velocity distributions of 
nucleons become isotropic in the central cell. Pions, however, kinetically equilibrate much later, at 
t = 8 fm/c. This effect is caused by the non-zero formation time for non-leading particles. Full local 
thermal equilibrium (LTE) in the central cell (i.e. consistency of the particle spectra and yields with 
T and ps extracted from {e, pb, Ps}) is first reached at t = 10 fm/c. To distinguish the later fully 
equilibrated phase from the earlier reaction stage in which only nucleons show kinetic equilibrium, the 
statistical model fits to the early reaction stage are denoted by open symbols whereas the fits during 
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Figure 5: UrQMD prediction for hadron ratios in Pb+Pb collisions at midrapidity (full circles). The 
open squares denote a thermal model fit to the UrQMD calculation. The fit yields a temperature of 
T=140 MeV and a chemical potential of //# — 210 MeV. 



the LTE phase are shown in full symbols. During the phase of LTE energy densities of ~ 350 MeV/fm 3 
are never exceeded, i.e. a density above which the notion of separated hadrons would lose its meaning. 
However, the main thermodynamic characteristics of the cell, T and fiB change rapidly with time. 
This clearly demonstrates that a fireball type description of hadronic matter is inadequate. 

Let us now turn to the hadronic final state of the heavy-ion reaction. Figure |S| shows a UrQMD 
prediction for hadron ratios in central Pb+Pb reactions at 160 GeV/nucleon around mid-rapidity (full 
circles). The UrQMD prediction has been fitted with a statistical model, yielding a temperature of 
T=140 MeV and a chemical potential of /ib = 210 MeV. However, this fit has been applied to the final 
hadron ratios after freeze-out. Here, the underlying assupmtion of the statistical model - namely a 
state of (global) equilibrium - is not anymore valid, since the break-up of the system and its freeze-out 
is governed by differences in the interaction properties (i.e. cross sections) of the individual hadron 
species. 

To study the breakup of the system in greater detail, let us turn to freeze-out distributions for 
individual hadron species: Figure ^ shows the freeze-out time distribution for pions, kaons, antikaons 
and hyperons at mid-rapidity in central Pb+Pb reactions at 160 GeV/nucleon. The distributions 
have been normalized in order to compare the shapes and not the absolute values. In contrast to 
the situation at 2 GeV/nucleon, where each meson species exhibits distinctly different freeze-out time 
distributions JDJ, all meson species here show surprisingly similar freeze-out behavior - the freeze- 
out time distributions all closely resemble each other. Only the hyperons show an entirely different 
freeze-out behavior - the situation is even more extreme in the case of the multi-strange Q, which 
exhibits a very sharp freeze-out time distribution, distinctly different from all other hadron species 
p3[ . Whereas the common freeze-out characteristics of the mesons seem to hint at least at a partial 
thermalization, the hyperons show that even at SPS energies there exists no common global freeze-out 
for all hadron species. The same observation applies also to the distribution of transverse freeze-out 
radii. Since these distributions have a large width, the average freeze-out radius clearly does not 
define a freeze-out volume and therefore estimates of the reaction volume or energy density based on 




Figure 6: Normalized freeze-out time distribution for pions, kaons, antikaons and hyperons. As with 
the freeze-out radii, the times for the meson species are very similar. The hyperons again show a 
different behavior. 



average freeze-out radii have to be regarded with great scepticism. The large width of the freeze-out 
distributions is supported experimentally by HBT source analysis which indicate the emitting pion 
source to be "transparent" , emitting pions from everywhere rather than from a thin surface layer ]30| . 

Figure displays the respective sources from which negatively charged pions freeze-out. Only 
inelastic processes have been taken into account. Approximately 80% of the pions stem from resonance 
decays, only about 20% originate from direct production via string fragmentation. Elastic meson-meson 
or meson-baryon scattering adds a background of 20% to those numbers, i.e. 20% of all pions scatter 
elastically after their last inelastic interaction before freeze-out. The decay contribution is dominanted 
by the p, uo and k* meson-resonances and the A1232 baryon-resonance - no weak decays have been taken 
into account in this analysis. However, more than 25% of the decay-pions originate from a multitude 
of different meson- and baryon-resonance states, some of which are shown on the l.h.s. of figure [7]; e. 
g. the two contributions marked p* stem from the P1435 and the P1700, respectively. 

The analysis of the pion sources is of great importance for the understanding of the reaction 
dynamics and for the interpretation of HBT correlation analysis results. The 20% contribution of pions 
originating from string fragmentation is clearly non-thermal, since string excitation is only prevalent 
in the most violent, early reaction stages. 

In summary, we have studied the evolution of relativistic Pb+Pb reactions at CERN/SPS energies 
from the early non-equilibrium phase through a stage of local thermal equilibration (in the central 
reaction cell) up to its final hadronic freeze-out. The importance of "partonic" degrees of freedom in 
the early reaction stage does not imply that an equilibrated Quark-Gluon-Plasma has been formed. 
In the UrQMD approach the "partonic" phase has been modeled as an incoherent superposition of 
non-interacting partonic constituents. Furthermore, these "partons" retain their original correlation 
into hadrons - deconfinement is not implemented into the present UrQMD approach. The leading 
(di) quark interactions (among each other and with fully formed hadrons) constitute an interacting 
"mixed phase" (for the constituent parton dynamics in this model, see, however [33, |27j). In contrast, 
parton cascades |29|] allow for interactions among the partons only, while hadronic final state 
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Figure 7: Pion sources in central Pb+Pb collisions at CERN energies: 80% of the final pions stem from 
resonance decays and 20% from direct production via string fragmentation. Decay-pions predominantly 
are emitted from the p and uj mesons and the A 12 32 resonance. 

interactions are to a large extent neglected. 

In the intermediate reaction phase, matter in the central cell can be viewed as hadrochemically 
equilibrated and exhibits an isentropic expansion. However, this equilibrium stage is limited only to 
the central reaction cell and breaks up in the late, dilute reaction phase close to freeze-out. The freeze- 
out of the system, which is governed by the individual hadron properties, has been studied in terms 
of freeze-out radii for different hadron species and a source analysis for the contributions of different 
microscopic processes to the final pion yield. A complex freeze-out scenario emerges with species- and 
momentum dependent broad freeze-out radius and time distributions. 
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